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and 
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University of Southern California 


SUMMARY 

This report presents an investigation of the problem of maximum likeli- 
hood detection of one of M Poisson processes in a background of additive 
Poisson noise. When the observables correspond to counts of emitted photo- 
electrons, the problem models a discrete version of a coherent M-ary optical 
communication system using photon counters in the presence of background 
radiation. Consideration is given to an average distance and a detection 
probability criterion. The advantages of an M-ary-pulsed intensity set 
(Poisson intensities wholly concentrated in a single counting interval) are 
demonstrated. The performance of such intensity sets is exhibited in terms 
of error probabilities, pulse widths, signal- to-noise ratio, and channel 
capacity. Behavior as a function of number of intensities M is also dis- 
cussed. By appropriate conversion these latter results may be used for de- 
termining power requirements in an optical pulse position modulation system. 


I. INTRODUCTION 

The application of detection theory to the field of optical communications 
has been a subject of increasing interest. Since the output of a photodetect- 
ing surface is often modeled as sequences of electron ’’counts 11 , and since 
optical photoelectrons have been readily accepted as obeying Poisson statis- 
tics, the analysis problem is basically one of signal detection involving 
Poisson processes. The problem was first formulated in this context by 
Reiffen and Sherman [ref. 1] , and further contributions were made by 
Abends [ref. 2] , Kailath [ref. 3] , and Helstrom [ref. 4] . In this report 
the general problem of M-ary detection based upon observations of events 
described by a time-discrete Poisson process is investigated. 

Though the formulation of the problem is of a general nature, the 
principal application is to optical communications, and the practical limits 
of such a system will govern the mathematical assumptions imposed. 
Consideration is given to the divergence criterion for detection and to a 
criterion of maximization of probability of detection, both readily accepted 
as suitable design objectives. The intensity set yielding optimal perform- 
ance based on special cases of these criteria is shown to be a special type 
of orthogonal intensity set, composed of M disjoint intensities wholly con- 
centrated in a counting interval. Previously, the superiority of this type of 
signal set in binary detection had been shown by Abends [ref. 2] , based on 
a signal-to-noise ratio criterion, and by Kailath [ref. 3], based on distance 
criteria. This report represents an extension of these results to M-ary 
Poisson detection. 


The formulation of the problem follows that of Reiffen and Sherman 
[ref. l] . The occurrence of events over an observed interval AT is said to 
obey a Poisson process if the probability of exactly k (an integer) occurring 
is given by 


P(k) = e' nAT • (1) 

The parameter n is the average rate of occurrence and is called the intensity 
of the process. The average number of events occurring is then nAT and is 
often called the level of the process. If the events occur over a sequence of 
AT intervals in which the intensity may vary from one interval to the next, 
but is constant over each interval, we have a discrete time-varying Poisson 
process. In photodetection each event corresponds to the emission of an 
electron, which occurs upon arrival of a photon, each photon having a fixed 
energy. The level is therefore proportional to the average energy received 
per interval, while the intensity n is proportional to the average power (see 
section V). Thus, constraints upon level and intensity in Poisson processes 
are equivalent to energy and power constraints on the incident radiation. 

In optical PCM communications, information is transmitted, as shown 
in Figure 1-a, by sending an optical signal intensity modulated with one of a 
set of possible intensities. The modulated signal is corrupted by fixed inten- 
sity background radiation during reception, resulting in a process, the inten- 
sity of which is the sum of each intensity. The output of the photodetecting 
surface at the receiver is then a time-varying Poisson process of electron 
counts having the received intensity. In an M-ary system the transmitter 
selects one of a set of M intensities for the optical process. The receiver, 
after photodetection, counts the number of electrons in each of M intervals 
A T seconds long and attempts to maximum likelihood detect which of M 
intensities is controlling the observed process. It is assumed that AT is 
suitably shorter than the inverse bandwidth of the intensities so that the inten- 
sity remains approximately constant over AT. In addition we assume the 
counting interval is exactly known at the receiver by a perfect synchronization 
link. Thus, the above system can be modeled by the block diagram in Fig- 
ure 1 -b. The input signal corresponds to a discrete Poisson process, while 
the interference appears as additive Poisson noise, (Recall that the sum of 
independent Poisson processes is itself a Poisson process having an intensity 
equal to the sum of the intensities. ) The assumption of additive Poisson noise implies 
that the optical bandwidth is in general wider than the receiver bandwidth. This is 
referred to as "classical” counting. To obtain deviations from "classical” counting 
takes a considerable effort and requires very special circumstances not normally en- 
countered. The model is idealized since other sources of interference, such as thermal 
noise and dark currents, are neglected. With this model the M-ary Poisson detection 
problem can be formulated as follows. Let a sequence of events obeying a discrete 
Poisson process occur over a sequence of M disjoint intervals A T where M AT = T 
and the count over each interval is independent of all others. Let the observed process 
be controlled by one of M possible intensity vectors n + n , q= 1 , 2 , . . . M, where 
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Figure 1. -- Optimum receiver 



The non-negative nqi is thus the intensity of n during the i interval. Under 
a fixed energy constraint for each signal, it is^ required that 


M 


Vn .AT 

Lu T 1 


N, for all q . 


( 3 ) 


The intensity vector n represents constant intensity background noise super- 
imposed upon the desTred intensity. Let the corresponding number of events 
occurring in the i^ 1 interval be k^. The problem then is to determine which 
of the possible intensity vectors n q is controlling the received Poisson proc- 
ess by observing the sequence of independent counts k = k^, k^, . . . 

Under a maximum likelihood detection criterion and a priori equilikely 
intensities, it is well known [ref. l] that the optimal test is to form the 
likelihood functions: 
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where 





Q qi 

!og J 

(n . + n ) AT] 

qi ° J 

(5) 


and s elect n 4 

n as the true intensity if no 

other likelihood function ex- 
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ceeds Aq(kj? If a°likelihood draw occurs (more than one A (k) is maximum), 
it is known that any randomized choice amongst the maxima^ is equivalent. 

In the following a purely random selection in the case of likelihood draws is 
used. Equation (4) can be interpreted as a crosscorrelation of k with the 
a .’s, an operation readily performed by a digital crosscorrelator [ref. 1] . 


4 


II. THE DIVERGENCE OF THE DETECTION TEST 


The divergence between two intensities n. and n of the above test is 


defined as 


-J 


-q 


v' 


where 


D. - E, /.(A. ) - E, t (A. ) 

jq k/j iq k /q n 


( 6 ) 


A. = A.(k) - A (k) 

jq J ~ q “ 

and E /.( A) is the conditional average of A with respect to k given the inten- 
sity nj/ J Abends [ref. 2] has shown that for M = 2 (binary detection) and the 
condition n- = 0, the divergence normalized by the variance of A is maximized 
by a "pulsea" type of intensity, where the level of the process is wholly con- 
centrated in a single- count interval. Kialath [ref. 3] has extended this result 
by showing that, under a total energy constraint, other suitable forms of 
"distance" are maximized by similar "pulsed 11 intensities. We extend these 
notions here to the M-ary case and the equal energy constraint of Eq. (3). 

The average divergence of an M-ary intensity set|n^J will be defined as 


D 


M 


i-ZI 


D. 


jq 


(7) 


Since E^^(k^) = (n^ + n^) AT, the average divergence becomes 


D = 


A T 
M 2 


Z I I <*ii- v{ lo8 bi + v iT ] 


j q 1 


log pn q . + n o )AT]j = Z M Z V° g [‘ + 


- Z Z n ji log 

j t q 


/ n .AT \ 


( 8 ) 


where K = n AT. The non-negativeness of the n.. and n allows one to 
o b ji o 

write: 
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D < 


2ATW 

M 

i j 



n. .log 
J 1 


1 + 


n.. A T 

-J1 

K 


2 

M 



n..AT 
_Ji 

K 



n.. AT 
J 1 


= 2N log (1 + |r) 




( 9 ) 


as an upper bound under the constrain of Eq. (3). However, the first equality 
holds if the second term in Eq. (8) is zero, requiring n.., j / q to be zero 
for all i at which n . is non- zero; that is, the intensities must be mutually 
disjoint. The second equality in Eq. (9) holds if n_ = N for one i and n.^ = 0 
for all other i. Thus, the upper bound for D occurs if the intensities of the 
set are disjoint and wholly concentrated in single counting interval. This is 
satisfied with the set 

N q - jN/AT 8. q = 1, 2, ...M (10) 

where 8 . is the Kronecher delta. The above represents an M-ary "pulsed" 
intensity ^et with each intensity occupying one of M intervals. It is signifi- 
cant to note that any disjoint intensity set, no matter how many intervals 
are used, yields the bound of the first inequality of Eq. (9), but only the 
pulsed intensity set of Eq. (10) yields the second bound. Thus, of all dis- 
joint intensity sets, only the pulsed set maximizes D, which immediately 
implies only M intervals are required for maximization with M intensities. 
Lastly, it may be noted that with an average energy constraint over all inten- 
sities: 


M Z X n ji AT = N 

j i 


instead of Eq. 


( 3 ), 


n.. < MN/A T 


results and Eq. 


(9) becomes 


(H) 


D < 2N log [l + ^~] (12) 

which exceeds that previously derived. Furthermore, the upper bound in 
Eq. (12) occurs when M - 1 intensities are zero everywhere, and one inten- 
sity is a pulsed intensity having a value MN/ AT. 
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III. DETECTION PROBABILITY 


The optimality of the M-ary pulsed intensity set has been shown, based 
on a divergence criterion. This section will show that in certain cases this 
superiority also extends over a criterion based on maximization of the detec- 
tion probability. First an expression for the detection probability for a 
general intensity set { n } is required. Usually this is obtained by first 
v writing the conditional pS-obability density of Aq(k), and then integrating over 
regions of correct decisioning. However, Aq (k) in Eq. (4) is a weighted sum 
of independent Poisson variates which in general is not a Poisson variable. 
Rather, the true density involves an M-fold convolution of modified Poisson 
densities, yielding a result that is difficult to integrate. Instead an alterna- 
tive expression for the detection probability is used. It was derived in the 
appendix and has the form: 

p d = ttZ < 13 > 

r m <1 

M 

where N is the intensity energy constraint of Eq. (3), R is the space of all 
M-dimensional vectors having non-negative integer components, and 

M f(n . + n )AT*1 

SMq.j) = n ~ . q . 1 , — ±_ exp f-Mn AT] . (14) 

i=l V L ° J 

The derivation of Eq. (13) follows an analogous procedure used in Gaussian 
channels (see Viterbi {ref. 5} ) , but is somewhat complicated by the fact that 
likelihood draws occur with non- zero probability. 

Determination of the intensity set{ n j for which P^ is maximum is of 
interest. This has been obtained for two^particular cases of interest: 

Case I — M = 2 and symmetric intensity sets: - 

Let M = 2 and consider the set of all possible symmetric intensity sets; 
i. e. , if nq = {a,b}, then n 2 = {b, a} . For this case it is easy tc^ show that 
for any intensity set of this type, the vector s_j for which 4' (l,j_) < (2,jJ 

when a>b is simply the set j = {jj, j^} such that jj£j 2 - Using the constraint 
of Eq. (3) and letting n^ = |a, N-aJ. and = |N-a,aj., -g- < a < N, the de- 

tection probability is: 
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where A = (a + K) and B = (N - a + K). Since A and B are positive, the 
above substantiates P D as a monotone increasing function of a. Therefore Pp 
is maximized with "a" having its maximum value, a = N, corresponding to 
the pulsed intensity set of (10) with M = 2. 
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Case II -- Any M, N/K — ► 0 


The limit above implies a high background noise level situation. We ob- 
serve here that 


- N 

D * V L m “ * <q ’i' 


-N 


R 


M 


M 


I 


R 


M 


C max exp 

q 


M f n .AT' 


i= 1 


> (16) 




N „ 

K“^° 


-N 

> 

M /-, 


M 

C max exp 

M q I A 

1 = 1 


j. n .AT/K 
i qi f 


where C = e FI (l/jjl ) and the limit follows since n .AT/K< N/K— 0. 

i qi 

Now, with the constraint of Eq. (3): 

V j.n .AT/K <j N/K (17) 

/ , J i qi * — J max 7 ' 7 


where 


i = max 

J max 

i 


Thus: 




-N 

P < — — 
D - M 


Z c {w N / K }' 


N 

K 


■=7- — - 0 


The upper bound occurs when 


max 

q 


j“P Z h “ qi * T / K j = «P {w N / K } 


(18) 


which clearly is true for the pulsed intensity set of Eq. (10), signifying 
asymptotic optimality for any M. 
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To determine optimal intensity sets (either global or local) in the general 
case, using Eq. (13), still remains an elusive task. It has been conjectured 
by many (e. g. , see refs, [l] and [3] ) that the pulsed intensity set is, in fact, 
the true optimal set, but to the authors’ knowledge a rigorous proof has not 
been shown. 


jf 
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IV. ERROR PROBABILITIES WITH PULSED INTENSITY SETS 




In this section the performance of the pulsed intensity set in M-ary de- 
tection is investigated by evaluating the error probability Pg = 1 - Pp>. This 
can be obtained by using Eq. (13), but the computation can be more con- 
veniently handled by noting that for the pulsed intensity set of Eq. (10) {A } 
of Eq. (4) constitutes a set of independent, Poisson random variables. Tne 
variable A q has level (N + K) if the q^ intensity is sent, and has level K 
otherwise (K = nAT). Recall, if the q^ 1 intensity is sent, a correct decision 
will be made with probability l/r + 1 if A q equals r other A' s and exceeds the 
remaining M - 1 - r. Therefore, when all possibilities are considered, the 
conditional detection probability is 


D/q 


-(N + MK) M- 


M 


x-1 

V K k - 
L TT e 
t=0 ‘ 


K 


^ ^ (N + K) X c -(N + K) 

L L L x: 


M - 1 - r 


-n r=0 x=l 


r (M - 1) i 

[ rj (M - 1 - r) J (r + l)j . 


-i r 


K 

x! 


-K 


(19) 


The right side is independent of q and thus represents the average detection 
probability. By applying the identity: 


M- 1 

I 

r =0 


(M - 1)! 

(r +1)1 (M - 1 - r)! 


M- 


Ur B r =££■ 


one can rewrite the error probability as 


P e (N, K, M) = 1 - p d = 1 - 


-(N+MK) 

Si 


00 p 

1 

x=l L 


(N + K) X e " (N+K) 


x-1 


y 

i—i f i 

. t=0 


t -K 


M-l 


1 

Ma 


(1 + a) M - l] 


( 20 ) 


where 


a 


K 



r 

4-J 

W 

& 
i 


x! 

Lt=o J 
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The parameter P p (N, K, M) has been digitally computed by the authors for 
various values of N, K, and M. An exemplary plot is shown in Figure 2 in 
which P E (N, 3, M) has been plotted for various M as a function of N. It is 
of paramount significance to note that P^ depends on both signal energy N 
and noise energy K in the counting interval, and not simply on their ratio. 
This fact is emphasized in Figure 3 in which P^ (N,K, 2) is plotted as a 
function of K for two fixed ratios N/K. This dependence on both signal and 
noise energies distinguishes the Poisson detection problem from the anal- 
ogous coherent gaussian channel problem. Note that the interfering noise 
energy K depends only upon the background energy in the interval A T, which 
is the width of the transmitted intensity pulse. The prime advantage of op- 
tical systems is precisely their ability to remove the effect of background 
noise by making A T small, and has been emphasized in the literature 
(ref. 6). * This fact can be illustrated graphically, using Eq. 20, by con- 
sidering a binary Poisson channel (M=2) sending information at a rate l/T 
bits/sec. The effect of the parameter A T is indicated by plotting 
P (N, n TA T/T, 2) as a function of A T/T, for fixed intensity energy N and 
background noise energy per bit interval n T. This is shown in Figure 4. 
The results indicate the continuous improvement obtained by decreasing the 
’ 1 duty cycle’ 1 AT/T, the ultimate limit corresponding to A T = 0. The 
improvement, of course, is made at the expense of information bandwidth 
and peak power (both inversely proportional to AT). Surprisingly, the 
improvement is quite small at low values of N, and the increase in band- 
width may not be worth the improvement obtained in error probability. 

A quantity of particular interest to communication engineers is the 
detected M signal-to-noise ratio”. This is often denoted [ ref. 7] as the 
ratio of the square of average electron count with no noise to the variance 
of the count when noise is present. For Ppisson counting statistics with 
pulsed intensities this becomes (S/N) = N / (N+K). The behavior of P^ of 
Eq. (20) as a function of K for fixed (S/N) is illustrated in Figure 5 for a 
binary system. The results again indicate the ambiguity in using S/N as a 
design criterion. The asymptotes show the wide functional variation of P E 
as K increases from zero. 

As illustrated in Figure 2, the error probabilities increase as M in- 
creases. However, the use of a single set of curves to compare various 
M-ary systems is misleading. An M-ary system with A T second counting 
intervals transmits log^M bits of information in MA T seconds. It therefore 
communicates at a rate 

io g m 

R = mIt bit/sec. (21) 


*Also Karp, S. , and Gagliardi, R. : A Low Duty Cycle Optical Communi- 
cation System. Presented at Eastcon Convention, Washington, D. C. , 
October 1967. 
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Figure 4. -- Error probability versus normalized pulse width 





If the transmission rate is normalized for each M, AT must be readjusted 
to maintain a fixed rate R = R . The effective noise level per counting 
interval is then 

/n \ 21og ? M 

, ”o AT = n o log 2 M/MR o = (j |J -Jjf- - 2K o log 2 M/M (22) 

where K is the noise energy in an interval l/2R . Thus, for a comparison 
of different M systems, each with fixed information rates, one should com- 
pare the parameter P_, (N, 2K log M/M, M) for each N. If this adjustment 
is made using Eq. (20y, the curves of Figure 6 are generated, with K q = 1. 

The curve corresponding to M = « is also shown, and is determined by 
taking the limit of (N, 2K log M/M, M) as M-* °o. This can be obtained 
by replacing K in Eq. (20) by K 1 = 2K o log M/M and noting: 

lim e -(N+MK'U 0 
M — oo 


lim 
M — 


(1 + a) M - 1 


oo 


Ma 


— * 1 for x > 1 


lim 

M -*■ oo 


fx-1 

I 

_t=o 


(K *) 4 -K' 


,n m - l 


0 for x = 1 

1 for x > 1 


Using the above results in 


lim 

M — oo 


P„(N, 2K log M/M, M) — 1 
hj o 


= 1 - (1 - Ne' N - e" N ) □(!+ N)e 


^ (N) X e -N 

L *' 

x=2 

-N 


(23) 


which is plotted as M = oo in Figure 6. It is noteworthy that Eq. (23) is 
precisely the probability of an event count of zero or one occurring in a 
noiseless counting interval of signal energy N. This has the following inter- 
pretation. As M -*■ oo the number of intervals becomes infinite, but the 
normalized noise energy per interval, K 1 = 2K log M/M, approaches zero. 
The probability that more than one event will occur in any one of M - 1 
independent non- signalling intervals having noise energy K' is given by 


1 


[1 + K'e 



M-l 
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Figure 6. -- Normalized error probability versus M 
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This approaches zero as M — », indicating that counts of zero or one will 
occur in every such interval with probability one. Furthermore, there will 
be an infinite nupn.be r of intervals with a zero count and with a one count. 
Therefore, as M-"**, an error will occur (with probability approaching one) 
whenever the signalling interval has a count of zero or one, and an error 
will never occur when the latter interval has a count greater than one; hence, 
Eq. (23). 


It is also interesting to note in Figure 6 that best system operation, in 
^terms of minimal error probability, does not always correspond to M-* <x > . 
In fact, it can be shown that best M operation depends strongly on the amount 
of background noise K . For example, if K ^ 0, it is easy to show X using 
Equation (20)] that for°M finite Pg(N, K', M) » (M - l)e /M, which is 
monotonically decreasing with M and always less than the M = °o value of 
Eq„ (23). Thus, with negligible background noise, system operation im- 
p^oves vfrith decreasing M and is best for M = 2. Physically, this means 
the noise reduction advantages due to decreasing A T as M increases does 
not offset the increasing errors due to the larger numbers of likelihood 
draws that will occur. (Recall a random choice is made in the event of 
draws. ) For large amounts of background noise, however, the converse is 
true, and M = <# does yield minimal error probability. 


It should be emphasized that a fixed energy constraint was imposed on 
the signal intensity, and therefore the time average power P q = N/T = 
NR/log M actually approaches zero as M “•’<» . If the average power level - 
of the source had been fixed at some level P , then N must be replaced by 
P log M/R in the previous equations, and it^an be seen that Pg -*■ 0 as 
hi -*• oo for any P > 0. This result may be compared to a similar - re suIF 
for an additive gsPussian channel (ref. 8) in which zero error probability 
occurred only if P satisfied a condition dependent on the rate R. 


The Pg results above are useful for determining the channel capacity 
(maximum information rate) of an M-ary pulsed intensity set. Assume a 
transmitter sends one of a group of M-pulsed intensities every T second, 
with each pulse having width A T = T/M. If the transmitter operates at a 
fixed rate R , then again T = log M/R as given by Eq. (21). The channel 
can now be. represented as a symmetric channel in which/ each of-the M 
equal likely intensities is converted to itself with probability 1 - Pg.^and is 
converted to any of the other intensities with ; equal probability 
Pg/ M- 1 . The channel capacity for this type of system is known to be 


C = 


log M + Pglog + (1 - P E ) log (1 


log M/R 


P E> 


(24) 


where Pg = P (N, n log M/MR , M). Again the signal intensity energy N, 
the backgroundnoise 0 power n , °and the rate R to be held fixed are con- 
sidered. Then, as M -*■ oo,Eg approaches the^imit in Eq. (23), while the 
channel capacity has the limil: 


r 


XT 


-N 


1 o 



for N finite. The above indicates information transfer can be forced to 
approach any desired rate with a finite signal energy by using an increasingly 
larger number of intensities and adjusting R at the transmitter. However, 
each level is transmitted with a non-zero error probability, and the infor- 
mation bandwidth and peak power become infinite. Again introduction of a 
transmitter power constraint, instead of an energy constraint, will yield 
operation at a capacity R q with a zero error probability, as M - *■ °o . 
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Y. CONCLUSIONS 


This report has investigated the M-ary Poisson detection problem, 
defined as maximum likelihood detection of one of a set of M discrete 
Poisson processes in the presence of an additive discrete Poisson noise 
process. The problem model represents a discrete version of an optical 
communication system in which the observables represent counts of 
photoelectrons, the signals represents intensity modulated CW optical 
sources, and the noise represents background radiation received within the 
V optical bandwidth. The photoelectron count can then be modeled as a time- 
varying Poisson process, the average rate of which is proportional to the 
sum of the intensities of the modulated source and the background radiation. 
In practical operation the intensity of the optical source is a continuous 
process, but the analyses may be put on a discrete basis by partitioning 
the signalling intervals into subintervals over which the intensity is taken 
to be constant. The above Poisson model is examined, and the advantages 
of a pulsed type of intensity set is demonstrated. The latter corresponds 
to an optical system using pulse position modulation in which information 
is transmitted by a burst or pulse of optical energy located in one of a set 
of pulse positions. The performance of such a system, in terms of pulse 
width and numbers of pulse positions, is presented herein. The results of 
this report basically represent theoretical limits to which an optical link 
can approach, since the deleterious effects of receiver (thermal) noise 
have been neglected. This latter assumption becomes valid, for example, 
when photomultipliers are used in detection, and the background radiation 
collected at the receiver is the predominant source of noise. 

The analyses and performance results are in terms of N and K, the 
average electron counts due to signal and noise, respectively. However, 
these results may be easily converted to average power requirements using 
the relations 

N = n P g M/hfB 

K = 77 P n /hfB 


where P and P are the average signal and background noise power, h = 
6.62 x 10-34, n Tj is the photodetector efficiency (including photomultipli- 
cation), f is the optical frequency of the CW source, and B = l/AT. The 
average power P and P can be further converted to transmitted power by 
introducing space losses n and receiver optics (e. g. , ref. 9, Ch. 1,2). Exact 
synchronization has been assumed here between transmitter and receiver at 
all times. In addition to receiver thermal noise, the analysis has excluded 
the effects of photomultiplier statistics, saturation, and dark currents. 
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APPENDIX 


This appendix presents the derivation of Eq. (13). The average proba- 
bility of correctly determining the true intensity in M-ary transmission is 

1 M 

P D = M £ P < D /“1> < A -‘> 

q=i 

where P(D/q) is the probability of correct detection, given that n + n is 
the true intensity. Now the conditional probability of the occurrence oi an 
observed vector k = j_ = {jj. j 2 » J 3* * - * * S iven 5. + n Q , is 


P(k = j/q) 


M 

n 

i=l 


[ (n qi +n o )A T ^i -(Mn )A T -N 
■ a -. ■ , e o e 

V 


(A -2) 


= ^ (q, 1) e" N 

where N is the energy constraint given in Eq. (3). The conditional detection 
probability P(D/q) is then obtained by summing over the set of all j_ such ^ 
that a correct decision is made. A correct decision will occur, when the q 
intensity is used, if Aqis selected as being the largest. If no other A^ ex- 
ceeds Aq but r of the A t 's equals Aq ? the receiver will be correct with a 
probability of l/r + 1, assuming a purely random selection is made when 
likelihood equalities occur. Nowj_is an M-dimensional vector with non- 
negative integer components, and the space of all such vectors is denoted 
as R^. The conditional detection probability P(D/q) can therefore be 
written by summing over allj_€ rM leading to a correct decision. Thus, 


P(D/q) 


M 


z 

r=0 


1 

r+1 



*(q,j_) 


-N 

e 


(A- 3) 


where J is the set of j_e R such that no other A^. exceeds A and r other 
A 's eq3al A . If we let I denote the r dimensional index set Corresponding 
tolhese r A^'s*, for simplicity, J can be denoted symbolically as 


J = / j CR^: A = max A 

q r q k 


A t' t 0 q} • 


(A-4) 
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Substituting Eq. (A-3) into Eq. (A-l) yields a general expression for the 
detection probability: 


-N M M-l 


— L ^ 1Y1 1Y1 — A 1 

P D = IT Y, X X r+T * (q ' ^ 
q=l r=0 J 


(A-5) 


Now by examining carefully the set J , the expression above can be simpli- 
fied, Making use of the monotonicitv q of the exoonential function, one can 


J qr = {l ^ rM; ex P ( A q ) = exp (max A fc ) = exp A t> t Cl q } 


TV K J-r j. > atj. • 

i CR : n [ (n . + n )A T] J i = max II [k.+ n )AT] ] i 

i = l qi o k i=1 ki o 


ni(n ..+ n )AT] J i, t Cl 

i=l q 


= CR M : ^ (q, j_) = max ^ (k,|) = ^ (t, j), t C I q j . (A-6) 

Thus Jq r can be alternatively defined as the set of j for which ^ (q, j) is one 
of r + l maximum (k, j) functions. This means every j_ in Jq r also belongs 
to r other sets Jtr»^CIq* or correspondingly, a point J_ in Jtr» tCIq, ^exists 
such that (t, jj = (q, j_). Note that the set of subspaces {Jqr} are dis- 

joint for different r but not for different q. With these facts consider the 
summation, 


Y* ^ * (q.j.) 

Z, Z, r+ 1 

q=l J 

A m r* 


(A-7) 


for fixed r. For any term of the sum, say SS r (q , j )/r + l, there exists r 
other terms having the same value, one for eac?i point j_ Q of , t Cl . The 
total contribution to the sum above from this set of r 4- l°terms r is theSi° 


" ®(q >3) I 

V [ — r ~ + I" J = ^ (q o’j 


o ) = max V (q, j Q ) , 


(A-8) 


the last equation following since j_ Q e Thus, overlapping points in the 

summation of Eq. (A-7) contribute a total amount given by Eq. (A-3). It 

therefore follows that 
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(A-9) 


M 

I I 

q=l J 


^(q.l) - 


r + 1 


qr 


= y 

yj max ^ (q, j_) 

„ 9 r q 


where O J is the union over of the subsets {j } . Inverting the order of 
summa*tioi£in Eq. (A-5) and using Eq. (A-9) allSws one to rewrite: 


V 


-N M-l 



-N 

^rr- 2 L max (q, jj {A- 10) 

r M q 


where we have employed the fact that |U Jqrj are disjoint subspaces, and 

the sum over all r spans the whole space R^. Equation (A-10) is Eq. (13) 
of the report. 
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